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We present a Molecular Dynamics study of the collective dynamics of a model for the fragile 
glass-former orthoterphenyl. In this model, introduced by Mossa, Di Leonardo, Ruocco and Sam- 
poli [S. Mossa, R. Di Leonardo, G. Ruocco and M. Sampoli, Phys. Rev. E 62, 612 (2000)], the 
intramolecular interaction among the three rigid phenyl rings is described by a set of force constants 
whose value has been fixed in order to obtain a realistic isolated molecule spectrum. The interac- 
tion between different molecules is described by a Lennard Jones site-site potential. We study the 
behavior of the coherent scattering functions Ft(q,t), considering the density fluctuations of both 
molecular and phenyl rings centers of mass; moreover we directly simulate the neutron scattering 
spectra taking into account both the contributions due to carbon and hydrogens atoms. We compare 
our results with the main predictions of the Mode Coupling Theory and with the available coherent 
neutron scattering experimental data. 

PACS number(s): 64.70.Pf, 71.15.Pd, 61.25.Em, 61.20.-p 



I. INTRODUCTION 

The understanding of the supercooled liquid and glassy 
phases in molecular systems is, nowadays, one of the main 
tasks of the physics of disordered materials (see and 
references therein for a general review). Two main the- 
ories provide us a description of the glass transition re- 
spectively from a thermodynamical and dynamical point 
of view. 

The first one (see Q and references therein) is based 
on first-principles computation of the equilibrium ther- 
modynamics of glasses and consider the glass transition 
as a true thermodynamic transition. In this context, the 
incoming of the glassy state is associated with an entropy 
crisis, i. e. the vanishing of the configurational entropy 
of the thcrmodynamically relevant states. 

The second approach is the Mode Coupling Theory 
(MCT) which studies the long-time structural dy- 
namics and its relation with the glass transition; in this 
context this transition has to be considered not as a regu- 
lar thermodynamical phase transition involving singular- 
ities of some observables, but as a kinetically induced 
transition from an ergodic to a non-ergodic behavior. 
The structural dynamics becomes so "slow" that the sys- 
tem of interest appeared frozen on the experimental time 
scales. 

On the experimental side, the collective dynamics of 
a huge variety of molecular systems has been investi- 
gated by means of several experimental techniques; col- 
loids |7|,|] and ortho-terphenyl (OTP) [g-[l2| are the most 
widely studied fragile Q supercooled liquids. They are 
only examples of an enormous experimental work (the 
interested reader is referred to Ref. 01 for an accurate 



comprehensive review). Moreover, in the last ten years 
the analysis of experimental results has been flanked by 
extensive use of numerical techniques, mainly Molecu- 
lar Dynamics (MD) and Monte Carlo (MC) simulations. 
The almost exponential growth of computational capa- 
bilities allows us to reach simulation times of the order 
of microseconds for simple monoatomics systems and of 
several tens nanosecond in the case of complex molecu- 
lar systems; such performances definitely permit com- 
parisons of numerical results with the structural very 
long time properties of the real systems. We remember, 
among many others, the results concerning the structural 
dynamics of two molecular liquids like SPC/E water Jp| 
and OTP (TJJT^]. We emphasize here the fact that the 
models involved in these studies are molecular and rigid 
in the sense that they have a structure and take into ac- 
count the orientational properties of the molecules but 
they disregard the role played by internal degrees of free- 
dom on the overall dynamical behavior of the system. 

In a recent paper Q we have introduced a new model 
for the intramolecular dynamics of orthoterphenyl (OTP) 
(T m = 329 K, T c ~ 290 K, T g = 243 K), one of the most 
deeply studied substances in the liquid, supercooled and 
glassy state. The introduction of such a flexible model 
allows us to understand the role of internal degrees of 
freedom in the short time {fast) dynamics ]l6| taking 
place on the time scale of a few picoseconds and to study 
their possible coupling with the long time center of mass 
dynamics ]T^ |. Moreover would allow us to emphasize 
once more the universality of the Mode Coupling Theory 
(MCT) approach for supercooled liquids and, in particu- 
lar, of its molecular version in the case of complex molec- 
ular liquids. In Ref. jl| we introduced in details the in- 
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termolecular model and we showed some results based 
on MD calculations; the self (one particle) dynamical 
properties of a system composed of 108 molecules have 
been studied in details and we have found good agree- 
ment with the main predictions of the MCT and with 
the experimental results related to the self intermediate 
scattering function ||l8|,[l9) and the self diffusion p0| , pl[ . 
In the present paper we complete the picture consider- 
ing the collective highly cooperative structural dynam- 
ics controlling the rearrangement of big portions of the 
system. We describe the dynamics of the molecules at 
the level of molecular centers of mass and at the level 
of phenyl rings center of mass. Moreover we calculate 
the neutron coherent scattering function taking into ac- 
count both the contributions due to carbon and hydrogen 
atoms in order to make a direct comparison with the ex- 
perimental results. 

The paper is organized as follows: in Sec. || we sum- 
marize the main predictions of the Mode Coupling The- 
ory (MCT) and we give a schematic introduction to the 
Molecular Mode Coupling Theory (MMCT); in Sec. |n 
we briefly describe the model and we recall some com- 
putational details. In Sec. IV we study the temperature 



and momentum dependence of the collective dynamics of 
the molecular and of the phenyl rings centers of mass. In 
Sec. we make a deep comparison between the exper- 
imental results of neutron scattering and the simulated 
neutron spectra calculated taking into account the inter- 
actions with both carbon and hydrogen atoms. Finally, 



Sec. VI contains a discussion of the results obtained and 



some conclusions. 



II. MODE-COUPLING AND MOLECULAR 
MODE-COUPLING THEORY 

As we have already stressed in Sec. | two main theo- 
ries, one intrinsically thcrmodynamical |4[ and the other 
purely dynamical |||(| have been developed up to now to 
describe the phenomenology observed in the glass transi- 
tion. In this paper we will make a comparison of our nu- 
merical results with the main predictions of MCT about 
the center of mass structural relaxation dynamics. The 
reason for this is twofold: first of all this is, indeed, the 
theory taking into account the states actually accessi- 
ble by the system on the time scale of a typical sim- 
ulation. Moreover, even if some experimental results, 
like the presence of the so-called knee characterizing the 
low-frequency behavior of the light scattering susceptibil- 
ity |2^,|23| or the presence of a cusp in the non-ergodicity 
parameter p4| , seem to contradict some of its predic- 
tions, such theory has been verified to hold for different 
experimental data on a very wide time region. Finally, 
in the last two years successful efforts have been made in 
order to generalize the theory to liquids of rigid molecules 
of arbitrary shape [E5J (see also Ref. p(| for the particu- 



lar case of a liquid of linear rigid molecules) taking into 
account both translational and rotational dynamics. 

In Ref. § we have summarized the main predictions 
of the so called ideal MCT so that here we recall only the 
fundamental equations of the theory, concerning the col- 
lective intermediate scattering function and that we will 
use in the following sections. MCT interprets the glass 
formation as a dynamical transition from an ergodic to 
a non-ergodic behavior at a cross-over temperature T c ; 
the theory is written as a self-consistent dynamical treat- 
ment of the intermediate scattering function , i. e. the 
time correlation of the density fluctuations of momentum 

This theoretical scheme can be considered the mathe- 
matical description of the physical picture of the cage ef- 
fect. Following the dynamics of a tagged particle it is pos- 
sible to recognize two main dynamical regions. On small 
time scale of order of some picosecond, (the /3-region) the 
dynamics of the particle is confined into a limited region 
(the cage) builded up by the nearest neighbors. In this 
regime it is possible to write the intermediate scattering 
function as 



F(Q,t) = f(q)+h(q)^¥-^G±(t/T <r ) (1) 

where f{q) is the non-ergodicity parameter (also referred 
to as the Debye- Waller factor), h(q) is an amplitude in- 
dependent of temperature and time and the ± in G± 
corresponds to time larger or smaller with respect to t ct , 
a parameter that fix the time scale of the /^-process. At 
this stage, the time dependence of the correlation func- 
tions is all embedded in the q-independent function G±. 
G±(t) is asymptotically expressed by two power laws, re- 
spectively the critical decay 

G + {t/T a ) = {t/T a y a T «t« T<J (2) 

and the von Schweidler law, 

G- (t/r a ) = -{t/T tT ) b r a « t « r a ; (3) 

characterized by the temperature and momentum inde- 
pendent exponents a and b; here r Q is the structural re- 
laxation time. 

At time scales of order of r a the cages start to break 
down and the particle starts to diffuse approaching pure 
brownian motion. This long time part of the dynamics is 
the so called a-region and is well described by a stretched 
exponential function 



F(q,t)~f(q) exp 



t 



(4) 



which verifies the time-temperature superposition princi- 
ple (TTSP). The a time scale r Q depends on temperature 
trough a power law of the form 
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r a oc (T - T c )~ 



(5) 
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The momentum dependence of the dynamical param- 
eters in the collective dynamics case is not trivial as in 
the single-particle case. When looking at the structural 
collective dynamics studying different values of the mo- 
mentum q means studying the highly cooperative time 
evolution of cages of average dimension 2it/q; it is clear 
that such time evolution is strongly coupled to the static 
topological structure of the system. More precisely, MCT 
predicts that the parameters f(q), h^ 1 (q), f3 a and r a (q) 
oscillate in phase with the static structure factor S(q). 

The result concerning the momentum dependence of 
the collective relaxation time is quite general and is well 
known as de Gennes narrowing |27j| . A general rela- 
tion |2£| holds among the one-particle t s and the collec- 
tive t c relaxation times, namely r c {q) ~ S(q)T s (q); if the 
diffusion limit is appropriate for F s (q,t), i.e. for values 
of q close to the first peak of the static structure factor, 
we obtain 



T c (q,T) 



1 5(g) 
D(T) q 2 



(6) 



where D(T) is the diffusion coefficient at temperature T 
and S(q) is supposed to be nearly temperature indepen- 
dent. 

In summary, the momentum dependence of the collec- 
tive dynamics is non trivially driven by the static struc- 
ture of the system. In particular, for molecular systems, 
the small length scale structure is determined by orienta- 
tional properties of the single molecules; a pure molecular 
translational dynamics will, obviously, loose all the dy- 
namical features controlled by the high momentum part 
of the static structure factor. 

An important step toward a correct explanation of 
the dynamics of molecular systems is to write down a 
Molecular Mode Coupling Theory (MMCT) H@ tak- 
ing into account both translational and rotational de- 
grees of freedom. If we consider a system of N identi- 
cal rigid molecules of arbitrary shape described by their 
center of mass positions fj (t) and by the Euler an- 
gles f2j(t) = {<pj(t),0j(t),Xj(t)) we can write the time- 
dependent microscopic one-particle density as 



N 



P (f, n,t) = J2 s[f - f n (t)]s[a, n n (t)}. 



(7) 



Expanding with respect to the complete set of func- 
tions given by the plane waves and the Wigner matrices 
D mn {Vl), we have the tensorial density modes 



Plr, 



M t) = J (21 + l)^ ^ fn{t) Dtn [«*(*)] • (8) 



Then, the generalization of the intermediate scattering 
function to the molecular case is the tensorial quantity 



These correlators are directly related to experimental 
quantities |26|; for I — V — 0, they describe the dy- 
namics of translational degrees of freedom which can be 
measured by neutron scattering when looking at the cen- 
ter of mass low frequency part of the spectrum; if the 
molecules possess a permanent dipolar moment, the cor- 
relators with I = I' = 1 give information related to 
dielectric measurements and I = I' = 2 is finally re- 
lated to the orientational contribution to light scatter- 
ing. At this stage, provided the static angular correlators 
Simn;i'm'ri {<}, 0) and the number density, it is possible to 
give a closed set of equations for the matrix S that com- 
pletely solve the problem of a liquid of rigid molecules. 

The problem is that MMCT seems to be not enough 
for a high structured molecular system like OTP; we will 
see that it is not possible to explain some features of the 
momentum dependence of the structural dynamics with- 
out taking into account the internal degrees of freedom 
(i. e. rotations of the side rings with respect to the cen- 
tral one) that turn out to be strongly coupled to the long 
time behavior of the density fluctuations. 



III. MODEL AND COMPUTATIONAL DETAILS 



In this section we give a brief description of the model 
and we address the reader to Ref. Q for details. 

In our model, the OTP molecule is constituted by three 
rigid hexagonal rings of side L a = 0.139 nm represent- 
ing the phenyl rings; two adjacent vertices of the central 
ring are bonded to one vertex of the two lateral rings by 
bonds of equilibrium length Lb — 0.15 nm. In the iso- 
lated molecule equilibrium position, the two lateral rings 
lie in planes that form an angle of about 54° with respect 
to the central ring's plane. In the model the two lateral 
rings are free to rotate along the molecular bonds, to 
stretch along the bonds and to tilt out of the plane iden- 
tified by the central ring. The intramolecular potential is 
then written as a sum of harmonic and anharmonic terms 
each one controlling one of these features. Every term is 
multiplied by a coupling constant whose actual value is 
determined in order to have a realistic isolated molecule 
vibrational spectrum. The intermolecular interaction is 
of the site-site Lennard- Jones type; each of these sites 
corresponds to a vertex of a hexagon and is occupied by 
a fictious atom of mass Mch = 13 amu representing a 
carbon-hydrogen pair. The actual values of the param- 
eters olj and clj have been fixed in order to have the 
first maximum of the static structure factor S(q) in the 
experimentally determined position p9| and to obtain 
the correct diffusional properties |20 2 
been fixed to the value r r — 1.6 im?. 



IF; the cut-off has 
It is worth noting 
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here that obviously the parameterization of the potential 
cannot be perfect; in our case it is possible to reproduce 
quite well the experimental results on the whole inves- 
tigated temperature range shifting the MD thermody- 
namical points at temperatures of 20 K above their true 
values. The MD simulated system is composed of 108 
molecules (324 phenyl rings for a total of 1944 Lennard 
Jones interaction sites) ; at each time step the intramolec- 
ular and intermolecular interaction forces are calculated 
and the equation of motion for the rings are solved for 
the translational and rotational part separately. 

Wide temperature and momentum ranges have been 
investigated for values of temperature 380 < T < 44QK 
and momentum 2 < q < 30 nm" 1 (the runs details are 
shown in Table III of Ref. (!]]) and the total simulation 
time is of almost one hundred nanoseconds. 



IV. MOLECULES AND PHENYL RINGS 

The collective density fluctuations dynamics is embed- 
ded in the coherent intermediate scattering function in 
general defined as 



F t (q,t) 



1 



NS(q) 



N N 

EE' 

i=l j=l 



-iq-[S,(t)-2 3 (0)] 



(10) 



where N is the number of molecules involved and S(q) the 
static structure factor. In the present case the position 
variables Xk{t) can be identified with different quantities; 
here we are interested in the dynamics of the molecular 
centers of mass and of the phenyl rings so that we will 
consider the following scattering functions 

F f {M) (q t) = - (V e -<9-[%(*)-Jfi- e «(0)]\ (ll) 



F t (R) ( q ,t) 



1 



N R SW(q) 



(EE 

ij c'S" 



-iq- R u ,(t)-R u „(0 



(12) 



Here M^> (t) is the position of the center of mass of the 
molecule £' at time t (£' = 1 . . . Nm), Ri£' (t) the position 
of the center of mass of the phenyl ring i (i = 1 ... 3) 
pertaining to the molecule £'; the functions are renormal- 
ized to the corresponding static structure factors. From 
now on, the superscripts (R) and (M) will be referred to 
molecular and rings quantities respectively. 

As in the case of the incoherent scattering function, 
at every temperature investigated we have reconstructed 
the whole curve, even on very short time scales, by mean 
of two sets of system configuration campionated with dif- 
ferent frequencies (see Table III of Ref. Q). At every 
investigated temperature, we considered the momentum 
values qi = 14 nm" 1 and (72 = 19 nm" 1 corresponding 
to the first and second peak of the experimental static 
structure factor, averaging on the values of q falling in 



the interval q ± Aq with Aq — 0.2 nm" 1 . Moreover, the 
momentum dependence of the principal dynamical pa- 
rameters has been investigated at T = 280, 300, 330 K 
for values of momenta ranging from 2 to 30 nm" . 

We made a long time analysis in term of the usual 
stretched exponential form Eq. (Q) determining the tem- 
perature and momentum dependence of the fitting pa- 
rameters r Q , (3 a and f q and verifying the TTSP. 

In Fig. |l] we show F^ R \q,t) calculated at qi for the 
temperatures T = 280, 300, 320, 350, 370, 390, 410, 430 K 
(from top to bottom); as in the case of the self dynamics, 
every curve decays to zero in the considered time window 
and it is clearly visible the two steps decaying pattern. 
The long time part of these F^R) have been fitted to 
Eq. (|4|) and the parameters r a , f3 a and f q are determined 
by a least square fitting routine. In the inset we plot the 
same curves as a function of the rescaled time t — tjT a ; 
all the curves collapse pretty well on a single master curve 
as predicted by the TTSP. The temperature dependence 
of the non-ergodicity parameter f q R ^ (top panel) and of 
the stretching parameter (bottom panel) are shown 
in Fig. |^; they are temperature independent, as predicted 
by MCT. The mean values f [R) = 0.78 and (3 {R] = 0.83 
(dashed line) have to be compared with the values de- 
termined in the case of the self dynamics f q ~ 0.7 and 
(3 a ~ 0.8 0]. The two values of (3 a are equal in the limit 

( R) 

of the error bars; at variance, the value for f q in the 
collective case is greater than the value found in the one 
particle case. 

In Fig. |§| we plot the structural relaxation times at qi 
and <j2 for molecules (triangles and diamonds) and phenyl 
rings (circles and squares) in order to test if both Tq M ' 1 
and follow the same power law which is supposed to 
be momentum independent. Both sets of data have found 
to be consistent with a power law of the form Eq. (|^) with 
parameters T c ~ 268, 7 ~ 2.3; these values have to be 
compared with the results concerning the one particle 
dynamics T c = 276 ±7 K and 7 = 2.0 ±0.4 0. The inset 
show the data plotted as a function of T = T — T c , in 
order to stress the power law dependence, and rescaled 
by an arbitrary factor in order to maximize the overlap. 

We now consider the momentum dependence of the 
collective dynamics at few selected temperatures T — 
280,300,330 K, spanning the momentum region in the 
interval 2 4- 30 nm -1 . We test the long time dynamics in 
term of the stretched exponential function and we verify 
the MCT predictions on the von Schweidler time region, 
characterized by the power exponent b. 

In Fig. U (top panel) we show the stretched exponential 
fits (solid lines) to F^ at some selected values of q; they 
works pretty good at least for time values greater than 5 
picoseconds as is more clear looking at the bottom panel. 
Here we show the whole calculated curve (circles) up to 
t ~ 0.05 ps at q = 14 nm" 1 and T — 300 K together with 
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the stretched exponential fit (solid line). From the top 
panel of Fig. || is qualitatively clear that the relaxation 
time depends non trivially on the momentum values. As 
we reminded in Sect. 0, MCT predicts that in the col- 
lective case, at fixed temperature, the relaxation times 
oscillate in phase with the static structure factor, i. e. 
they are strongly coupled to the static structure of the 
system. 

In Fig. ^ and Fig. ^ we plot the momentum depen- 
dence of the collective relaxation times (top panels) for 
molecules and rings respectively, at T = 280, 300, 330 K 
(circles, squares and triangles respectively) together with 
the corresponding static structure factors divided by q 2 
(bottom panels). The correlation among these two quan- 
tities, predicted by Eq. (^|), is evident; it is also evident 
that, at variance with the temperature dependence, the 
momentum dependence of the relaxation times is com- 
pletely different in the two cases. 

In the molecular case (Fig. |J) only a maximum at 
q ~ 9 nm -1 corresponding to intermolecular correla- 
tions is present ( a small shoulder at q ~ 14 nm _1 , re- 
lated to correlations between rings pertaining to differ- 
ent molecules, can be also identified). In the case of the 
phenyl rings (Fig. |^) the momentum dependence is much 
more structured and three main features are present: 

(i) a maximum at q ~ 9 nm -1 related to correlations 
between molecular centers of mass 

(ii) a shoulder at q ~ 14 nm -1 (well developed in a 
maximum at the lowest temperature T = 280 K) re- 
lated to correlations between rings belonging to different 
molecules 

(iii) a maximum at q ~ 22 nm" 1 related to correlations 
between rings pertaining to the same molecules; this is 
the more important features related mainly to the ori- 
entation of the lateral rings with respect to the central 
ring. 

It is clear from this result that in complex molecu- 
lar glass formers there are intramolecular rotational and 
vibrational degrees of freedom that couple to the transla- 
tional long-time modes as already stressed in Ref . ; a 
theory not taking into account these degrees of freedom 
cannot explain the whole momentum dependence of the 
centers of mass dynamics. 

In Fig. [?] we show the different relaxation times (left 
panel for molecules and right panel for rings) at the 
three selected temperatures multiplied by the corre- 
sponding diffusion coefficients from Ref. |j| ; as predicted 
by Eq. (^), these products are temperature independent 
for values of q close to the first maximum of the static 
structure factor. 

In Fig. ^|(A) we plot the momentum dependence of the 
stretching parameter (3^ while in Fig. 0(B) we show the 

(R) 

non ergodicity parameter fg at the three selected tem- 
peratures together with S^ R \q)/q 2 ( Fig. ||(C)); also in 
this case a strong correlation between the oscillations of 



the quantities is clear, as expected. 

We have seen in Sec. || that the long time limit of the 
/3-region can be described by the von Schweidler power 
law Eqs. (0) (||). The exponent b is expected to be mo- 
mentum independent and to assume the same value of 
the self dynamics case namely b — 0.52 at the con- 
trary, the amplitude h(q) is expected to be momentum 
dependent and to oscillate out of phase with the static 
structure factor. We then calculated a power law fit in 
the form F t {R} [q, t) = f (R) - c (R) (q) t bW for all values of 
momentum considered, in a time window depending on 
the particular q- value but always included in the interval 
2 4- 30 ps; moreover we considered the three parameters 
free as in the case of the self dynamics. All the observa- 
tions done in the previous work concerning the great un- 
certainties on the estimated values of the fitting param- 
eters, hold in the present case. In Fig. [9] (A) we plot the 
power exponent b(q) that is supposed to be momentum 
independent; some smooth oscillations are nevertheless 
present but this can be due to the interplay during the 
fitting procedure with the other oscillating parameters. 
In Fig. |9|(B) we plot fq R ^; these points must be compared 
with the results of Fig. || (B) and the agreement looks 
quite nice as expected. The value of the plateau, indeed, 
must be the same if determined as the small time limit 
of the a-region or the long time limit of the /3-process. 
In Fig. H (C) we plot the quantity l/c^iq) and also in 
this case oscillations in phase with S R (q)/q 2 are evident, 
as expected. 

V. NEUTRON SCATTERING 

Neutron scattering is one of the most powerful tools 
used in the study of supercooled liquids and glasses in the 
q— region covered by the MD simulations. Experimen- 
tally the scattering function F t (q, t) of Eq. ( |l0| ) can be 
determined JTcj ] by neutron scattering experiments either 
directly, on neutron spin echo instruments, or Fourier 
transforming the dynamical structure factor S{q, uS) 

calculated by means of triple axis backscattering or time 
of flight spectroscopy. The experimental neutron scatter- 
ing cross section (da/dftdE) is generally composed of a 
coherent and an incoherent part 

da/dVtdE ~ < b > 2 S coh (q, w) (14) 
+ [< b 2 > - < b > 2 }S mcoh {q,uj) 

where b is the scattering length and the symbol () denotes 
an average over the distribution of nuclear spins and iso- 
topes. The isotopic composition of the sample allows 
us to study selectively the collective motion via coher- 
ent scattering from deuterated samples (the scattering 
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lengths of D and C atoms are basically coincident) and 
the one particle motion via incoherent scattering from 
protonated samples. 

The interaction of neutrons with a bulk sample of OTP 
can be simulated numerically taking into account the in- 
teractions of neutrons with both carbon (C) and deu- 
terium (D) atoms. H atoms are not considered in our 
dynamics but is a reasonable approximation to put them 
in fixed position on the line extending from the center of 
the ring trough a carbon atom at the fixed C-H distance 
dc-H = 0.107 nm; so that, knowing the coordinates of 
the rings, is trivial to reconstruct their own positions. 
We then define a neutron (N) coherent scattering func- 
tion Ff N \q, t) as 

w A' A" ij 

where Na = Nc + Nh (3456 in this case) is the total 
number of atoms, r\>i£> (t) is the position of the atom A' 
pertaining to the ring i in the molecule £' and S^ N ' (q) 
is the static structure factor of Fig. 11 of Ref. Q). The 
number of hydrogen atoms is 4 for each central ring and 
5 for each lateral ring. The scattering lengths 6a are in 
principle different for the carbon and deuterium atoms 
but, as we observed in Ref. [Q, they are both positive 
and of the same magnitude so that is a good approxima- 
tion to consider the product by b\n an ineffective positive 
constant. The function F} (q, t) is the quantity directly 
comparable with the experimental data. 

In the present section we present a comparison be- 
tween the temperature and momentum dependencies of 
the MD and experimental spectra of Ref. Jl2] calculated 
from perdeuterated C18-D14 by means of coherent neu- 
tron time-of-flight and backscattering spectroscopy. 

At this stage few remarks must be done about the mo- 
mentum dependence of the MD and experimental sets 
of data. All these data are supposed to depend on the 
structure of the systems so that, in general, some differ- 
ences are expected (experimental and MD structures are 
slightly different JjJ). Anyway, two observations about 
the experimental results must be made. First of all, as 
reported in Ref. p2[ , some reservation is necessary for the 
experimental data at the smallest momenta, q < 6 nm -1 , 
where f q tends toward 1. Indeed, in this region one ex- 
pects significant background from incoherent scattering, 
which contributes about 15% of the total cross section, 
and from multiple scattering. Moreover the technique 
used to determine the values of the dynamical parame- 
ters from the experimental data are quite different with 
respect the MD computation. Indeed, due to the lim- 
ited dynamical window of the available spectrometers, 
in the experimental case a direct fit of the data to the 
stretched exponential function with three independent 



parameters is not possible. Based on the observations 
that r(q) cx r](T)/T (r](T) is the viscosity at tempera- 
ture T) and that the line shape is independent on tem- 
perature, at fixed momentum q the spectra at different 
temperatures are rescaled in time to t — t/t s where the 
scaling time is given by t s = t v (T)/t v (T = 290K) and 
T n = T)(T)/T. In this way the data converge toward 
a temperature-independent long time asymptote; this is 
the curve actually fitted to the stretched exponential. 

In Fig. [h] we show the temperature dependence of the 
stretched exponential fitting parameters for the a-region 
of both MD (open circles) and experimental (solid circles) 
data determined at q = 14 nm -1 . In panel (A) we show 
the structural relaxation times; the inset shows the same 
data but the MD points have been shifted of 20 K as al- 
ready explained in Ref. Q. We stressed out, indeed, that 
a non perfect parameterization of the diffusive behavior 
of the model system controlled by the value of clj, has 
shifted the MD thermodynamical point about 20 K above 
the corresponding experimental temperature. From the 
figure is clear that the agreement among the two sets of 
data is very good on a very wide time region. In panel 
(B) we show our results for the non ergodicity parameter 
f q and also in this case the agreement is very good among 
the two sets of data; in panel (C) we plot the stretching 
parameters (3 a . In this case the MD points are system- 
atically above the experimentally determined data. This 
effect can be due to the observations made before; more- 
over, for the determination of the f3 a parameter, the very 
long time points are crucial and they seem to lack in the 
experimental data analysis. 

With Fig. [n] we start the comparison of the momen- 
tum dependence of the two sets of data. In the top 
panel of Fig. |ll| we plot the simulated scattering func- 
tions F^ N) (q,t) (symbols) at T = 290 K at the in- 
dicated q together with the stretched exponential long 
time fits (solid lines). In this case fq is expected 
to be momentum dependent so that, in order to have 
a master plot (bottom panel of Fig. [ll]), we have to 
show Ff(q,t) = Ff(q,t)/fq as a function of the 
rescaled time t — t/ri N \ The rescaled curves collapse 
quite well on a master curve as expected. In the next 
three figures we show the comparison for the momen- 
tum dependence of the MD data at the three tempera- 
tures T — 280, 290, 300 K and the experimental data at 
T = 313, 320, 330 K; the two sets of temperatures should 
be comparable due to the 20 K shift of the MD data. 

In Fig. [l2| we show the MD (left panel) and experi- 
mental (right panel) relaxation times together with the 
corresponding static structure factors renormalized to q 2 . 
Striking similarities are clear even if maxima on MD re- 
sults correspond to bumps of experimental results. The 
clearest difference is the decrease of t« at small value of 
q that appear to be completely absent in the experimental 
case. This could be due to the incoherent background on 
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the experimental data stressed above; moreover, is well 
known that sometimes MD is not able to determine the 
correct value of r at small value of q but usually it tends 
to overestimate its correct value, at variance of this case. 
Anyway, the agreement among the two sets of data at 
the higher values of q corresponding to the intramolecu- 
lar correlations is surprising; none rigid model could be 
able to do that. 

In Fig. [l^ we show the results concerning the non er- 
godicity parameter fq N \ The different curves are tem- 
perature independent as expected; both sets of data show 
a maximum at q ~ 14 nm -1 , while the maximum at q ~ 9 
nm -1 for the MD data correspond to a little shoulder in 
the experimental case. Also in this case a decreasing part 
at small q is present in the MD data at variance with the 
experimental results. But according to Fig. 7 of Ref. 
the calculation of the non-ergodicity parameter fitting 
the /3— region gives a plateau at q ~ 6 nm -1 at variance 
with the increasing behavior of the a— region analysis. 



However, as already stressed in Sec. |IV|, these data are 
supposed to coincide. Concluding, also in this case, the 
experimental data at small q seem to be not reliable. 

In Fig. [li] we finally show the stretching parameter 
/3a . MD data are very noisy, nevertheless it is possi- 
ble to recognize an oscillatory behavior; moreover, the 
points at the smallest available value of q seem to catch 
the decreasing behavior of the experimental data. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have concluded the analysis of the 
long time center of mass dynamics of the intramolecular 
model for OTP introduced in Ref. there we found 
good agreement among the diffusion properties of the 
simulated and the real system as well as among the two 
single-particle dynamics. Moreover in Ref. M we found a 
very good agreement with the main predictions of MCT. 
In the present paper we definitely confirm such an agree- 
ment. 

We have studied the collective density fluctuations on 
a large temperature and momentum range range, con- 
sidering both the fluctuations related to molecular and 
phenyl ring centers of mass. With respect to the tem- 
perature dependence we found the usual double step de- 
caying pattern and we confirmed the main predictions of 
the MCT about the behavior of the stretched exponen- 
tial parameters; in particular the relaxation times obey 
the same power law. The momentum dependence of the 
stretching parameters appears really interesting; MCT 
predicts that the behavior in the momentum space of 
the structural relaxation time, at variance with the the 
trivial square law of the one-particle case, is driven by 
the structure, namely it is proportional to S(q)/q 2 . In 
both cases, molecules and phenyl rings, this prediction is 



completely fulfilled. The phenyl rings behavior is partic- 
ularly interesting; every feature of the quantity S(q)/q 2 
is mirrored on the tA (q) curve at the three investigated 
temperatures T = 280,300,330 K up to a value of mo- 
mentum q ~ 30 nm . In particular, the maximum at 
q ~ 22 nm -1 is related to fluctuations taking place on 
molecular length scales; the long time structural dynam- 
ics appears, indeed, coupled to the dynamics of internal 
molecular degrees of freedom. Similar oscillations are 
present for the other stretched parameters and also the 
correct behavior for the short time /3-region is found. 

The next step has been a comparison among the ex- 
perimental and simulated neutron scattering spectra cal- 
culated considering both the scattering from carbon and 
deuterium atoms. The temperature dependences of the 
relaxation times in both cases are well described by the 
same power law for almost three decades in time; a good 
agreement is also found for the other stretching param- 
eters. The momentum dependence is qualitatively very 
similar in both cases; although some different features 
not clear are present. 

At this point some conclusions can be drawn about the 
capability of our model to describe the long time dynam- 
ics of the real system. All the center of mass time scales 
calculated by molecular dynamics have been found to be 
consistent with a power law, although some discrepan- 
cies were present for the actual value of the power expo- 
nent mainly due to uncertainties on the fitting procedure. 
Moreover it has been shown that our actual MD thermo- 
dynamical point is shifted by about 20 K with respect to 
the corresponding experimental point. Taking into ac- 
count all these informations, we plot in Fig. [l^ all the 
time scales related to the centers of mass dynamics con- 
sidered up to now, both experimental and numerical, as 
a function of the rescaled temperature T — T — T c where 
T c = 290 K H| for the experimental points and T c = 270 
K for the MD results. In particular, tv-f{T) (circles) is 
the time scale related to the shear viscosity r/ s of Ref. [ pO[ , 
r a (T) (left triangles) is the neutron scattering structural 
relaxation time of Ref. E3| at q = 14 nm -1 and the in- 
verse of the experimental diffusion coefficient (squares) 
is from Refs. (2(]] |^l| . The MD one-particle relaxation 
time calculated on phenyl rings tJ (T) (diamonds) and 
self-diffusion coefficient (up triangles) are from Ref. [@; 
finally, ri R \T) (down triangles) is the structural relax- 
ation time of Fig. || calculated at q = 14 nm -1 on phenyl 
rings and Tq^' (T) the simulated neutron scattering re- 
laxation time of Fig. |To|. All the data collapse pretty well 
for almost three decades in time on the same line corre- 
sponding, on a double-log scale, to a power law of the 
form of Eq. (||) of exponent 7 = 2.55 |Tgj| . 

Concluding, our model has found to be a very success- 
ful model for the centers of mass dynamics of the real 
system, showing a critical behavior consistent with the 
experimental results in a wide time window. The imple- 
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mentation of an intramolecular dynamics is relevant to 
such extent; in particular, the internal degrees of freedom 
appear to be strongly coupled to the long time structural 
dynamics as is also clear from the study of the rotational 
properties of the system |3lJ. The next step will be the 
understanding of the role played by the internal degrees 
of freedom in the fast relaxations observed experimen- 
tally @. 
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FIG. 1. Intermediate 
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coherent scattering 



functions 



Fj: (q,t) calculated on phenyl rings at q = 14 nm for 
the temperatures T = 280, 300, 320, 350, 370, 390, 410, 430 K 
(from top to bottom); in the inset we show the same curves 
rescaled as a function of t/ri R K 
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FIG. 2. Temperature dependence of the stretched ex- 
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non-ergodicity parameter fq R \T). Bottom: stretching pa- 
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FIG. 3. Temperature dependence of the structural relax- 
ation times at q\ — 14 nm -1 and §2 = 19 nm -1 calculated 
both on rings (circles and squares respectively) and molecules 
(triangles and diamonds respectively) centers of mass. In the 
inset the date are shown in a double-log scale as a function 
of the rescaled temperature (T — T c ); the points have been 
shifted in order to maximize the mutual overlap and to stress 
the power law behavior. The power law of exponent 7 = 2.3 
is also shown (dashed line); the value for T c is 268 K. 
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FIG. 4. Top: Intermediate coherent scattering functions 
F} R ^ (<?, t) calculated at fixed temperature T = 300 K for dif- 
ferent values of momentum q; the corresponding long-time 
stretched exponential fits are also shown (solid lines). Bot- 
tom: Here we show the whole time behavior of Fj; R \q,t) at 
q = 14 nm -1 together with the stretched exponential fit. 
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tively). Bottom: Structure factor calculated on molecular 
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